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ABSTRACT 

The links between the internal structure of galaxy groups and clusters and cosmolog- 
ical parameters are reviewed here. The mass density profiles of clusters, inferred from 
both optical analyses of the galaxy surface number density profile coupled with inter- 
nal kinematics, and from weak gravitational lensing, appear cuspy with inner density 
profiles at least as steep as r~^, as is expected from high-resolution cosmological TV- 
body simulations starting with scale-free or CDM initial conditions. The high level of 
substructure both in the galaxy distribution and in the diffuse hot gas seen in X-rays 
is consistent with dynamically young clusters and a cosmological density parameter 
f^o ?J 0.5. Despite the importance of cluster-cluster merging witnessed by these sub- 
structures, the spherical top-hat cosmology appears to give good indications of the 
underlying physics of clusters and groups, including a fundamental evolutionary track 
in a space of observational parameters for groups. Dynamical processes operating in 
groups and clusters are reviewed. Previrialization of groups and clusters before their 
full collapse is not expected on theoretical grounds, and the ensemble of positions of 
compact groups, relative to the fundamental track is a potentially useful constraint on 
previrialization. 

1 Introduction 

This review focuses on the links between the richest hierarchies, groups and clusters, 
and global cosmological properties of the Universe, principally the density parame- 
ter, Qq, and the primordial density fluctuation spectrum, P{k). We begin with the 
observational facts, derive constraints on cosmological parameters, and focus in the 
end on dynamical and cosmological evolution of structures. The review is written for 
non-experts, so the specialist should be patient with some very obvious definitions 
given here. We will not touch upon issues of the long term evolution of galaxies in 
groups/clusters or the baryon content of these systems. 

2 Observational Facts 

2.1 Small-scale structures of the Universe 

The attractive nature of gravity pushes matter to concentrate in overdense regions, 
up to large scales, because of gravity's long range. To first order, galaxies are a tracer 
of the total mass content of the Universe (departures of which are called bias). And 
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so, one expects and indeed sees that galaxies tend to agglomerate in a hierarchy of 
structures as shown in Table |T[ 



Table 1: Hierarchy of galaxy systems 



System 


iVgal 


"^faint 


Isolation 


Scale 


{M/L) 


p 

^ gal 


Clusters 


30-300 


mg 


+ 


2 


None 


1.5 h-^ 


Mpc 


300 h 


10% 


Loose Groups 


3-30 


nil 


+ 


3 


None 


1.5 


Mpc 


150 h 


50% 


Compact Groups 


4-8 


nil 


+ 


3 


3 r 


< 0.2 h- 


^Mpc 


50 h 


0.1% 


Binaries 


2 


nil 


+ 


3 


5ri2 


< 0.2 h- 


^Mpc 


50 h 


10% 


Isolated 


1 
















30% 



Notes: A^^gai is the number of galaxies per system within the given range of magnitudes 
(m = — 2.51ogio /^wx + est). M/L is the mass-to-light ratio, in solar units, assuming 
dynamical equilibrium, Pgai is the fraction of galaxies in a given type of system, 
and h = ifo/(100 kms~^ Mpc~^) is the dimensionless Hubble constant, measuring the 
current rate of expansion of the Universe. The criteria for clusters |]1|, compact groups 
[0], and binaries [0], have been slightly modified here. 



Originally, galaxy systems have been cataloged through eye-ball searches of galaxy 
positions projected on the sky. Spectroscopy of galaxies provide redshifts, measuring 
galaxy radial velocities (through the Doppler effect), which are the sum of the Hub- 
ble expansion velocity (proportional to distance through the Hubble law v = HqD) 
and the peculiar velocity relative to this Hubble flow. Except for nearby galaxies, 
the Hubble flow term dominates, so that redshift is a first-order distance estimator. 
Recent catalogs of galaxy systems are defined in 3D (projected positions and depth 
inferred from the redshift) in automated fashion with well defined criteria. 

One should understand that galaxy systems probably span a continuous range of 
hierarchies, and therefore not put too much weight in the special categories defined 
here, whose properties depend strongly on the adopted selection criteria. From Table 
1^, most galaxies lie within fairly rich environments, and we shall see in § ^]l] that 



they are dense enough to be, using cosmological terms, in the quasi-linear to non- 
linear regimes of the growth of the primordial density fiuctuations by gravitational 
instability. Their properties are thus function of both the underlying cosmological 
parameters and the dynamical processes operating within them. 



2.2 Density profiles 

Surprisingly little is known on the galaxy number or total mass density profiles in 
groups and clusters? The galaxy number density profiles of groups and clusters are 
poorly known because of 

• small number statistics (most clusters have at best a few hundred galaxies with 
radial velocities confirming their membership) 
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• their unknown center 

• contamination from foreground and background interlopers 

For example, in clusters, model profiles with homogeneous cores fit the data 
but so do models with cuspy cores 0, and the two are hard to distinguish given the 
limits on the data |^. Homogeneous cores are the best fits when one superposes the 
information from different clusters, rescaling to the same size 0, but when one uses 
the X-ray barycenter or the position of the giant cD galaxy as the cluster center in- 
stead of the center of the galaxy positions, the clusters unequivocally show cuspy cores 
[^]. So homogeneous cores were the result of poor choices for the cluster centers. The 
cuspiness of clusters has been recently confirmed in an analysis of 1500 suspected 
members of a single nearby rich cluster (Coma), for which the number space-density 
profile is 1/ ~ r~^. The presence of possible interlopers in the sample can only wash 
out the central density cusp, so that, in reality, the slope of the 3D number density 
profile may be even steeper. 

If the situation in galaxy clusters is difficult, it should be impossible to conclude 
anything on groups of galaxies. Nevertheless, by averaging over groups in the NBG 
[|] catalog, and rescaling the groups to some unit size, it has been noted that 
these were closer to homogeneous than to ~ r~^. Averaging over the best defined 
compact group sample 0, again with rescaling, compact groups catalog were found 
to be significantly centrally concentrated [|T^, although the slope of the underlying 



profile is not given. A reanalysis of that compact group sample, after removal 



of groups with interlopers in redshift space |]12[, and starting with the distribution 
of absolute pair separations, shows that the distribution of galaxy separations within 
compact groups is consistent with a unique absolute density profile, falling off as 
in the envelope and with a small homogeneous core of size 18 kpc, that 



-2.4 



-P— = P Z2 ' 



is half the median pair separation within compact groups. It is not yet clear that 
this result is caused by high central concentration or by the presence of tight binaries 
within less concentrated density profiles (as seems to be the case in loose groups 0). 

What about the total mass density profiles? One can resort to either optical data 
or X-ray data. In both cases, one writes the equation of hydrostatic equilibrium (or 
Jeans equation), which in its general spherically symmetric form is 

dP^_d^^ GMtot(r) 
dr dr 

where p is the density (in mass or number) of whatever tracer we use to measure the 
pressure P. It is then easy to obtain the mass density profile 

1 rfMtot 

P = o — ; — ■ (^) 

4:7rr^ dr 

Now, from equation (|l|), the Jeans equation for the 'gas' of galaxies in a spherical 
cluster is 

2,g_!^_,M. (3) 

dr r r 
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where u is the space number density of galaxies in the cluster, cr^ is the radial velocity 
dispersion (standard deviation of the velocity distribution) of the cluster, and /5anis = 
1 — cr^/af. is the velocity anisotropy, given that at is one tangential component of the 
cluster velocity dispersion. 

The problem with the optical analysis is that the variation of anisotropy with 
radius is unknown, even if there are good reasons to believe that, in the central 
regions of clusters, the system should be isotropic, since the two-body relaxation 
time of galaxy-galaxy encounters is short (see § § below). The anisotropy of the outer 
envelope can be either radial, if the envelope evolution is dominated by near-radial 
infall of galaxies, or tangential if the cluster is experiencing an off-center collision with 
a smaller one. Note that the reversed problem of fixing the mass density profile and 
computing the anisotropy profile is solvable in quadrature |jl3[ . 



The advantage of the X-ray method is that there is no pressure anisotropy term 
to include in the Jeans equation. Using the equation of state of an ideal hot X-ray 
emitting gas, the total mass profile is in 



M (r) - ^"^^ fd\nn ^d\nT\ 
° Gfinip \d\nr d\nr J ' 

where n and T are the local gas number density and temperature, while firrip is the 
mean particle mass in the hot plasma. Since both the galaxy system and the gas 
reacts to the same overall potential, and when both systems are isothermal, and 
when the galaxy system is also isotropic, one finds [|15] (eqs. 0] and ||^) ngas ~ n'^si, 
where (3 = a1/{kT / fimp) is the ratio of kinetic energies of the galaxies and gas. Now, 
using equations (^ and (||), one can compare the total mass density profile to the 
gas density profile, where the emissivity of the gas locally scales as the square of its 
density. It turns out that the dark matter is significantly more centrally peaked than 
the diffuse X-ray emitting gas [|16|, |l^ |18|, |T^. However, the dark matter density 



profile is similar to the galaxy number density profile |]T^ |T^. 

The optical analysis of the Coma cluster, assuming isotropy of the velocity dis- 
tribution, yields [0 a total mass density p ~ steeper than the galaxy number 
density given above. Putting in radial anisotropy at large radii produces a broken 
total mass density profile varying as in the center but only in the envelope 0. 
It seems difficult to have the inner mass density profile fall as slowly as the galaxy 
number density profile with an anisotropy profile consistent with dynamical and cos- 
mological principles. If the Coma cluster is indeed regular, the X-rays may be telling 
us something about the radial variation of the velocity anisotropy. 

Gravitational optics is rapidly becoming a promising alternative method for de- 
riving projected mass distributions. Indeed, the amplification due to an elliptical 
gravitational lens is PU 

1 

^ ~ 7i ^2 2 ' 
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where the matter (or Ricci) term k, can be expressed in terms of a critical surface 
mass density T,c- 

1^ = ^ . (6) 

^2 1 



AnGDiil-Di/Ds) ' 

where Di and are the distances to the lens and to the source, respectively, and 
where the shear (or Weyl) term is 

^ = n- ^{x)xdx . (7) 

Equations (|^) and (^ indicate that gravitational amplification of background galaxies 
by the potential of a galaxy group cluster is directly related to the surface mass 
density. 

When the surface mass density of the cluster is close to Sc, background galaxies 
are lensed into a giant arc, arclet, or are just slightly elongated, in a sequence of 
decreasing alignment of the background galaxy with the caustics in the plane of the 



foreground cluster . While giant arcs have been observed P2| , P3| , there only comes 
one or two per cluster and thus cannot constrain the cluster density profile. On the 
other hand, the elongations of lensed background galaxies, caused by the shear term 
(eq. 0), allow one to map the surface density of a cluster, by measuring and inverting 
the shear field 0|. With such methods, one finds cuspy surface density profiles at 



least as steep as S ~ r ^ p5| , p6| . Surprisingly, there often seems to be more mass 



at a given radius than inferred from the X-rays (eq. by a factor 2 to 3 [^, ^ 



but the discrepancy disappears if one allows for small-scale fluctuations in the radial 



variation of the gas temperature |£6| 



2. 3 Substructure 

The naive view that clusters of galaxies are regular structures with little substructure 
can be taken as justification for the dynamical analyses of their underlying mass 
content. Recent studies point to the opposite picture of very irregular clusters. These 
are based upon analyses of the projected distribution of galaxies in clusters without 
[p9| , |30| , |31| or with |^ the radial velocities of the member galaxies, as well as on X- 



ray analyses. In particular, at least 30% of clusters show significant substructure 
from optical 3D analyses |3^, including the once canonical regular Coma cluster p5 |. 



Similarly, the X-ray analyses indicate that at least 22% of clusters show bimodal 
substructure [^], and optical analyses indicate that roughly half |Q or as much as 
90% 1^ of the clusters that are regular in the X-rays show small-scale substructure 
(the latter analysis is based upon wavelets). 
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In loose groups, again the small-number statistics make it very difficult to assess 
statistically significant substructure. In a system of say 8 galaxies, small-scale sub- 
structure would show up as galaxy pairs (subgroups of triplets or quartets of galaxies 
would show up as separate groups in the catalogs). For the NBG |^ group catalog, 
1.5 pairs are found on average per group [^], although a fraction three times lower 
has been claimed p8| for groups in an unpublished catalog. 



2.4 Fundamental plane 

Three parameters are most easily measured in astronomical systems: angular size, 
velocity dispersion (or temperature), and flux. In a subsample of systems at known 
distances, angular size provides the physical size, and flux provides the (intrinsic) 
luminosity {i.e., power). Thus, astronomical systems are well characterized by the 
parameters, R (size), V (velocity dispersion), and L (luminosity). Whether these 
systems are star clusters, galaxies, or galaxy clusters, people have attempted to find 
a plane in R, V, L space in which the objects of a sample all fit. Such fundamental 
planes have indeed been found for galaxy clusters ||3^ with the relation L ~ i?"'^ V^'^, 
not very different from that of elliptical galaxies (with no disks or spiral structure) 
or globular star clusters. 

Fundamental planes are useful for two reasons: 1) They teach us about how the 
systems form and evolve. 2) They serve as distance indicators, independently of 
the Hubble law (K = HqD, where Vr = cz is the radial velocity, while c is the 
velocity of light, and z the redshift). This is important, because the Hubble law 
is only uniform on average on very large scales, but on smaller scales is perturbed 
by local overdensities of matter (where the expansion rate is lower than the average 
one, and becomes even negative when matter falls into a galaxy, group, or cluster) 
and underdensities (where the expansion rate is larger). Therefore, outliers from the 
fundamental plane are usually systems for which the distance was poorly estimated, 
and the correct distance is obtained by forcing the system to lie on the plane. In 
other words, if L ~ CstR'^V'^ (no relation with the previous as and /3s), then the 
correct distance to an outlier is 

where the angular size 9 = R/D, and the measured flux / = L/^AnD^). 



2.5 Internal kinematics 

The motions of galaxies in clusters probe the potentials of these systems, when these 
are near dynamical equilibrium. For example, Jeans' equation (P) specifies the radial 
variation of the second moment of the velocity distribution, the squared velocity 
dispersion. The measured line-of-sight velocity dispersion is the emissivity-weighted 
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average of the velocity dispersion, where the geometry becomes a httle comphcated 
in the presence of velocity anisotropy |TB[. In any event, the radial variation of the 



line-of-sight velocity dispersion provides the knowledge of the velocity anisotropy, 
assuming that one knows the potential ||T3|, or of the potential if one knows the 



anisotropy. 

The early realization that the velocity dispersion decreases towards the center of 
the Coma cluster |^ (as well as outwards in the envelope) is interpreted as either 



a cuspy total mass density profile with a core that has isotropic velocities |^ or a 



homogeneous core with tangential (nearly circular) velocity dispersions A recent 
systematic study of a large number of clusters reveals that roughly 15% of clus- 
ters have such inverted velocity dispersion profiles, while roughly 40% have velocity 
dispersions decreasing outwards everywhere, the remaining profiles being either fiat 
or unclassifiable. The clusters with cuspy velocity dispersion profiles indicate cuspy 
steep density profiles (p ~ R'"^ , where 7 = 3.5 — 4) [Q, with at best very small 
(< 40 kpc) homogeneous cores . 



3 Cosmological Constraints 

3.1 The spherical top-hat approximation 

Consider a constant finite spherical overdensity in a homogeneous Universe. The 
evolution of this top-hat perturbation is governed by the simple equation 

ii = -GM/R^ . (8) 

It will expand with the Hubble expansion then reach maximum expansion at the 
turnaround time and collapse at twice the turnaround time. Subsequently, the mat- 
ter may reach dynamical equilibrium, and become virialized as it satisfies the virial 
theorem 2(T) + (W) = 0, where T and W are its kinetic and potential energies, 
respectively. 

It is easy to show (without cosmological arguments, simply integrating eq. [§) 
that the turnaround time can be expressed in terms of the density at turnaround 

rta = [-^) , (9) 

where pta is the density at turnaround. Once the top-hat collapses, the surrounding 
regions will collapse onto it in a gradual fashion, each reaching maximum expansion 
following equation (Rl) where pta is now the mean density within the turnaround region 



at the epoch of turnaround This process of so-called secondary infall has been 
studied in detail and for Q = 1, one obtains simple scaling laws for the matter 
within the turnaround radius, beyond which matter is expanding: 

Mt,~t^/3 i?ta~t^/^ Pt.-t-\ (10) 
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Reversing equation (^) one can then obtain the present radius of turnaround of 
cosmological objects in the top-hat approximation, and these are given in Table |^ 
from 

fSGMtX 



1/3 



R 



ta 



(11) 



independent of 43 . Here, to is the age of the Universe 



Table 2: Turnaround radii 



System 


Mta(M0) 


Rt. (Mpc) 


Small galaxy 


3 X 10^^ 


0.4 


Large galaxy 


3 X 10^2 


0.8 


Small group 


3 X 1013 


L7 


Small cluster 


3 X 10^^ 


3.6 


Large cluster 


3 X 10^5 


7.8 



Within the turnaround radius, matter is infalling up to the region of shell crossing 



inside which matter is well mixed (see for a graphical example). The radius of 



mixing is approximately the radius of second turnaround or rebound radius, which is 
roughly one-third of the turnaround radius For an Q = 1 universe, the density 



is critical at pc = {GnGt'^)'^ , and the mean overdensity within the turnaround radius 
is (using eq. [0) 97r^/16 = 5.55. If matter virializes at collapse and settles, from the 
virial theorem, at half its turnaround radius at epoch 2 Tta, its density will be 8 times 
larger than it was at turnaround, while the density of the Universe will decreased by a 
factor (eq. |T^) 4 (for Q = 1), so that the overdensity at virialization is 18 tt^ = 178. 
Loose groups of galaxies are defined with density contrasts of 20 |^ to 80 and so 



must still be feeling the collapse. Combining the virial theorem with this cosmological 
minimum density yields a small maximum size of ~ 200 kpc and minimum velocity 
dispersion of ~ 250 km s"^ for groups [0. 

Recent cosmological coUisionless A^-body simulations with power-law primor- 
dial density fluctuation spectra P{k) ~ k"' indicate that the spherical top-hat ap- 
proximation is very good, in that the radius where the density is 178 times the mean 
density delimits the infalling region from the quasi-static region. Nevertheless, the 
top-hat approximation cannot represent all the cosmological aspects of clusters and 
groups, because the Universe is filled with structures and substructures (§ |27^ ). 

3.2 Density profiles 

What should one expect for the density profiles of clusters? In simple terms, the 
final density profile of an object will be set in part by whether the object collapses at 
once or in a slower smoother fashion. The details are set by the primordial density 
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fluctuation spectrum, i.e., the variation of tlie rms overdensity fluctuations as with 
( CO moving) scale R. 

In the top hat approximation (§ it is easy to derive density profiles using the 
relations of equation (0), yielding p ~ r~^/^ and ~ [^|^ ^j^g —9/4 slope of 



the density profile has been verified with more detailed semi-analytical calculations 
|i5[] and with iV-body simulations . 

Consider now a spherical top-hat perturbation in an open (fi < 1) Universe. 
Because the Universe is less dense than critical, there will be a lack of matter falling 
onto the perturbation at late times, so that the envelope of the evolving system should 
be steeper than p ~ r~^/^ [^, as checked with iV-body simulations However, 



the regions that are accreted at early times came from a Universe whose density was 
closer to critical, thus resembling more an f2 = 1 Universe. 

Now reality is not as simple as spherical top-hat perturbations in an otherwise 
homogeneous Universe. First, perturbations may not be spherical, which then pro- 
duces shallower slopes [^. Moreover, fluctuations in the primordial density field 
occur everywhere on a variety of scales. One can then show analytically that the 
density profiles have slopes of Min[— 2, —3(3 + n)/{4 + n)], given a scale-free primor- 
dial density fluctuation spectrum P{k) ~ /c", and again, this has been checked with 



A^-body simulations |57, 54, 55]. A^-body simulations with CDM spectra (which have 
spectral index n increasing with scale from —3 to 1, and is roughly —2 on the scale of 
galaxies and —1 on the scale of clusters) show a varying a slope from —1 inside to —4 
or —3 p^l outside (the latter study includes dissipative gas). Interestingly, high- 



resolution scale-free simulations |^ reproduce the density profile with slope varying 



from —1 to —3 found in low resolution CDM simulations with gas All the high 
resolution simulations with CDM or scale-free spectra show cuspy cores, with a core 
radius smaller than the numerical resolution, which can be as small as a few kpc |l58| . 

3. 3 Substructure 

In a low density Universe, structure formation tends to freeze at redshift z = l/^o 
|60[| , because there is simply not enough matter to keep accreting onto already formed 
structures. In contrast, a critical density Universe sees accretion at all times. There- 
fore, the level of substructure in galaxy clusters seems to be an excellent way to 
constrain the density parameter Qq. 

The approach is to predict the fraction of galaxy clusters with given level of sub- 
structure, as a function of Q and compare with the observations (§ p. 31) . One way 
is the analytical computation of the distribution of collapse times of structures of 
given mass Alternatively, one builds Monte-Carlo realizations of the history 

of structures and their mergers, based upon the Press-Schechter mass function 



of structures at given epoch |4| (this is, in essence, a different extreme repre- 



sentation of the Universe than the top-hat model, which is limited to the accretion 
of numerous small objects onto groups and clusters). All studies yield Qq > 0.5, 
although the precise lower limit to the density parameter is still subject to debate. 
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These predictions are insensitive to the primordial density fluctuation spectrum |63 
and to the cosmological constant . 

For example, bimodal substructure is expected for a critical density universe in 
15% 24% and 28% depending on how long the bimodal substructure 



can survive within the cluster (assumed at 0.08 times the present age of the Universe, 
for the first two studies and 0.1 to, respectively), compared to > 22% observed (§ p.3| ). 
For f^o = 0.2 the corresponding predicted fractions are 4 to 6 times lower, and clearly 
discrepant with the observations. 



Half of all clusters are predicted to harbor 10% level small-scale substructure i6| 



if fio = 1, which compares well with the smaller of the two observed estimates |^ 



If instead 90% of all clusters harbor small-scale substructure [ETfl, one would need 



either a higher Qq or a survival time longer than 0.08 to- A recent study |^5|] does 
provide a 50% longer survival time for small-scale substructure caused by accreted 
small clusters, but slightly less, when the small-scale substructure if caused by the 
dense cores of two similar mass clusters detaching from their parents when these two 
merge together 



iV-body simulations with gas [^, ^ confirm the semi- analytical predictions given 



above and illustrate very well how the gaseous distribution in clusters in a low-f2 
Universe is considerably smoother on small-scales than observed in the X rays. 

3.4 Previrialization 

Most analytical approaches to cosmology apply the top-hat approach, in which the 
collapse is radial. In principle, non-radial motions can arise as angular momentum can 
be pumped into the system by external tidal torques or by gravitational interactions 
on small scales within the system. The existence of significant non-radial motions is 
known as Previrialization [^. How important is this effect? The first numerical study 



[ [fO|] argued for a very strong effect from the outer environment, so that the critical 
density for collapse, expressed in units of present day linearly evolved primordial 
density, is found to be of order 5, instead of the canonical value of 1.69 for fio = 1 Pll 



or 1.61 for ^o = 0.1 [Q. A more realistic set of cosmological A^-body simulations 
with P{k) ~ k~^, showed that power on small scales within clusters does not alter 
its collapse. Moreover, implicit in this study is that external tidal torques bring only 
low angular momentum (J = 0.07 Jcirc) just within the turnaround radius of clusters. 
A very recent analytical study [^], using perturbative expansions in the quasi-linear 



regime, suggests that for n ~ — 1, the collapse is as one expects from linear theory, 
hence with insignificant non-radial motions. The cold dark matter spectrum, still 



considered to be a close approximation to the true power spectrum |7J, has a slope 
close to —1 for the infalling regions of groups and clusters, hence previrialization 
should not be important in modeling these systems. 
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4 Dynamical Evolution 

There are quite a few dynamical processes competing in the evolution of groups and 
clusters of galaxies. What is written below is largely taken from a previous review 

The simplest timescale in a self-gravitating system is the circular orbital time: 

37r 



2nR 

K;irc(-R) 

since 1/^ = GM{R)/R. 



Gp{R) 



1/2 

(12) 



4-1 Relaxation 

As a test particle undergoes scattering collisions within a sea of field particles, it 
will progressively forget its initial conditions. This two-body relaxation time can be 
defined in at least three ways: 



1 dV'^X ' / 1 dE\ ' /rfsin^Aa^ ^ 



T"2-rel = 



dt 




where Aa is the deflection angle in an encounter. This can be written as |[76|] 

^3 

T2-rel 



G^mfUfiv/a^) In A 



where v is the velocity of the test particle, mf,n, and are the mass, number density, 
and ID velocity dispersion of the fleld particles, respectively, f{x) is the fraction of 
particles traveling faster than xa^, and In A, is called the Coulomb logarithm, with 
typical values of 2 to 10, where A is the ratio of maximum to minimum impact 
parameter. For a system of galaxies and dark matter particles, one flnds that the 
galaxies relax by galaxy-galaxy collisions, but not by collisions with individual dark 
matter particles (whose masses are too low). Similarly, the dark matter particles 
relax mainly by collisions with individual galaxies. 

A collective relaxation time has been derived, not by summing up the encounters 
but by computing the collective response of the system 



TAT-rel — Cst 



Grrifn 



2/3 



This collective relaxation turns out to be somewhat more efficient than two-body 
relaxation in clusters and loose groups but not in dense groups. In general, only the 
cores of rich clusters are relaxed. 
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Particles that evolve in a rapidly time varying potential can rapidly forget their 
initial conditions fTH] This violent relaxation occurs in a timescale 



,50, , 

Tv-rel ~ Tff ~ Tdyn whcn |-^| > |v ■ V0| , 

where Tg is the free-fall time, and is the global potential. This applies for example to 
collapsing systems, as is often the case in cosmology. For example, elliptical galaxies 
are thought to form by dissipationless collapse or by mergers of disk galaxies |^ 



and since both phenomena perturb the potential in a rapid violent manner, the cores 
of elliptical galaxies appear relaxed although their 2-body (and collective) relaxation 
times are much longer than the age of the Universe. 

Energy exchanges during encounters lead to energy equipartition, and the more 
massive objects tend to move slowly and fall to the center of their systems, leading 
to mass segregation. The timescale is similar to that of 2-body relaxation. 

Occasionally, collisions will pump sufficient energy into a particle that it's ve- 
locity will be larger than the escape velocity of its system, and, barring subsequent 
encounters that may reduce its kinetic energy, the particle will escape. If during 
one relaxation time the distribution of particle velocities reaches a Maxwellian, the 
fraction of unbound particles is roughly 1% so that the timescale for evaporation is 
roughly 100 times the relaxation time [pO|. 



4-2 Dynamical friction, orbital decay and circularization 

Field particles are deflected by the mass of the test particle. Hence, in the frame of 
the test particle, the field particle density is higher behind the test particle than in 
front of it. This leads to a drag force known as dynamical friction [^, which plays a 
major role in group and cluster dynamics. For an infinite homogeneous medium, the 
timescale for dynamical friction is 



_ / 1 cif||\ 

'^^'^\^\'dt) " 27rG2(m + m/)p/(t;/a,)ln(l + A2) ' ^^^^ 
where p is the local mass density of field particles, and where / and A are defined in 

§13 

Perhaps more physical is the timescale for orbital decay defined as 
(idRY^ (RdE/dR\ 3 /p l\ 

Dynamical friction and orbital decay can lead to misleading answers: 

• No orbital decay is predicted in zero density environments, whereas a satellite 
galaxy sitting just outside its parent galaxy will see its orbit decay, because of 
resonances with its parent 
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• Although orbital decay should be slowed by tidal effects that reduce the test 
particle's mass, the contrary may occur with a satellite galaxy circling its parent, 
as the tides from the latter remove stars from the former, and these carry off 
energy and angular momentum, thus accelerating the orbital decay 

• For a particle radially falling into a medium with a outwards decreasing density 
profile, the dynamical friction time computed from equation (|T3p is longer than 
the time on which the particle sees an increasing density. A system, dense 
enough to survive the tidal forces from the primary into which it falls, will 
nearly stick to the core of the primary after its second passage, if its mass is 
more than 10% of the primary's mass fS^ . 

For circular orbits, combining equations ([T^), (|TB|), and (P^, one obtains 

Tod ^ + Ip/p) M/m 2 M/m 

Torb 27r/(l) ln(l + A2)(l + m/M)V2 " vr (1 + m/M)V2 in [1 + (M/m)2/3] ' ^ ^ 

where / ~ 1/2, pmin ~ p^ax ~ R, primary and secondary have the same mean 
density {[M/m] ~ [R/r]^), and for singular isothermal density profiles [p ~ r'"^), 
p = 3p. The term (1 + m/My/'^ is a correction for non-negligible secondary to 
equations (|12D and (p!3|) . Table ^ below gives the ratios from equation (|l^) for typical 
astrophysical ratios. 



Table 3: Number of circular orbits for orbital decay 



Secondary 


M/m 


Torh / Tod 


Galaxy 


10 000 


1000 


Group 


100 


20 


Small cluster 


10 


3.5 


Moderate cluster 


3 


1.5 



Notes: Examples of secondaries are for a rich cluster primary, but should be consid- 
ered as academic, since the secondaries are not expected to acquire sufficient angular 
momentum as to orbit in near circular orbits as they enter the cluster ( §3.4| ). 

From table ^, dynamical friction and orbital decay are not important, in terms 
of number of near-circular orbits, unless the mass ratio of secondary to primary is 
^ 0.1, i.e., small clusters falling into larger ones, or galaxies falling into small groups. 
In absolute terms, the circular orbital time is r ~ 0.4to for a secondary orbiting a 
IO^^Mq cluster at 1 Mpc from its center. Hence, the timescale for orbital decay 
in rich clusters is greater than the age of the Universe for systems with m < 10^"^ Mq, 
orbiting at 1 Mpc. 
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Another outgrowth of dynamical friction is orbital circularization, whose timescale 
can be defined as the rate a test particle acquires angular momentum from interactions 
with other particles: 



1 dJ\ ^ 



Tn. 



Jarc{E) dt 



and is found to be shorter than the orbital decay time outside of the core of a 
cluster. 



4.3 Tides 



Tidal forces act on particles in a system relative to the full system itself. As such 
there are two types of tides acting on galaxies in groups and clusters: those caused by 
close encounters with other galaxies and those caused by variations in the gradient of 
the global group/cluster potential. The first type of tides {coUisional stripping) has 
a timescale 

/ 1 dm\ / N , w-i Cst 

Trs = I I = {[Am/m)n{(jv)) = — - — , 



m dt 



nr„v, 



where a is the collisional stripping cross-section, and the outer stars are assumed to 
follow elongated orbits ISf 



Global potential tides depend strongly on the galaxy's orbit around the cluster. 
If the galaxy is phase locked in a nearly circular orbit around the cluster, it will feel 
a roughly constant tidal shear, and its tidal radius will be obtained by equating the 
tidal shear at a given radius in the galaxy with the gravitational pull that the full 
galaxy exerts on a star at that radius, plus an inertial term: 



yielding for r <C -R 



(GM{r) 



Gm{r) 



(16) 



Pci(i?) , V^{R) 



+ 



;i7) 



i.e., the galaxy is tidally truncated at a radius where its mean density is of the 
order of the mean cluster density within the radius Rp of closest approach of the 
galaxy (where Vp and Kirc are the pericentric and circular velocities, respectively). 

If the orbits are elongated, the instantaneous tide obtained from equation ([T6|) is 
short lived and the galaxy experiences a tidal shock [^. Using the impulse approx- 
imation [89|, in which the perturber moves with a constant relative velocity V, one 
can show Q], again for r -C i?, that 



Pgin) = Cstpd(i?p)/(e) 
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where Rp is the pericentric of the galaxy's orbit, and /(e) is a function of order 
unity of the galaxy's orbital eccentricity. This criterion is similar to that for circular 
orbits, but the constants are higher, because at given pericenter, a galaxy in a circular 
orbit must feel a more effective tide, since it is long-lived |^]. Numerical simulations 



[pl[] confirm this result although other simulations suggest that the tide is most 
efficient for some intermediate elongation at given pericenter, when this is within 
the nearly homogeneous region of the cluster. Note that the timescales for global 
potential tides are basically the orbital timescales divided by the typical mass-loss 
per passage through the cluster core. 

The effectiveness of a tide is related to the maximum strength of the tide times 
the duration of this maximum tide. So, from equation ( ]T7| ) one gets 

A ^ A, - A. 2 - 3pcl/Pcl + ^//KL 

Av ~ FtidAt ~ pgAt - 




3 ( 1 - ^ 1 - ( 1 - 3^ I ( - 1 ) for > Kirc . 



Hence, the A^-body simulations showing an intermediate orbit elongation for maxi- 



mum global potential tides are understood, since when the cluster region is nearly 
homogeneous, the effective tide increases with increasing pericenter velocity, but not 
when the cluster density profile decreases sharply as outside the core of the Modified 
Hubble model used in the simulations. 



4-4 Mergers 



Because galaxies have their own internal energy, galaxy collisions are often inelastic 
enough to lead to merging. The timescale for merging may be estimated from a 
merging cross-section, again as 



Using a numerical experimental cross-section ||93|, and integrating over a Maxwellian 



velocity distribution, the merger time can be written pi 



Cst 



nrlvgK{vd/vg 



where n is the number density of galaxies, Vg and Vg are the galaxy half-mass radius 
and internal velocity dispersion, respectively, and the dimensionless merging efficiency 
K is optimum for groups {vd — Vg), while for clusters it falls off as v^i [0. In groups 
as dense as compact groups appear to be, merging ought to be extremely effi- 
cient, and the relatively low fraction of ellipticals indicates that chance alignments 
are contaminating the compact group catalogs |^4[. Despite their high velocity dis- 
persions, rich clusters are able to produce the right amount of mergers to produce 
elliptical morphologies [0], and moreover, merging is able to account [0] for the 



morphology- density |^| and morphology-radius |95| relations 
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5 Cluster evolution 

5.1 Dynamical evolution 

The physical processes described in § ^ compete in the evolution of the galaxy system. 
For example, merging leads to increased merger cross-sections, hence to a merging 
instability However, this instability is slowed down by tidal processes which are 
usually thought to truncate galaxies of their outlying particles which become unbound 
Yet, if the merging cross-section is related to the galaxy half-mass radius ||97|| , 



and since the tidal processes for galaxies on elongated orbits or from collisions pump 
energy into the system, then the half-mass radius of those particles that remain bound 
to the galaxy should increase. The question remains whether the new half-mass radius 
is then greater or smaller than the old value, but this reviewer is not aware of any 
numerical study that has addressed this question yet. 

Similarly, if galaxies possess huge halos when they enter clusters (as one can infer 
from Table |^), they should feel strong dynamical friction and orbital decay. But once 
they pass near the cluster core, these halos should be severely stripped by the global 
tide of the cluster potential, after which the galaxy is less massive and will be no 
longer subject to much dynamical friction. 

The interplay of the various dynamical processes and the difficulty in analytical 
modeling of tidal effects render necessary to run numerical A^-body simulations to 
see how groups and clusters evolve. In coUisionless particle simulations, galaxies 



rapidly dissolve in clusters a problem known as overmerging , which is often 
ascribed to numerical effects. In particular, two-body relaxation between galaxy 
particles and cluster particles has been blamed for overmerging Two analytical 



estimates differ on the importance of this effect. Whereas, one [ |10CI|| finds a short 
timescale for the two-body relaxation between galaxy particles and cluster particles: 
rci_g = (fg/fci)^Tci-ci — 0.1 Tci_ch the other study ||101|| considers the time r = U/U = 



[//(A^Af/imp.apx.), where U is the internal energy of the galaxy and N is the rate of 
collisions, and finds a timescale typically longer than the Hubble time. This last study 
[|101|| also checks that evaporation is not responsible for overmerging. 

Overmerging can also be ascribed to physical effects such as the tides from the 
global potential of the cluster or dense group or two-body relaxation between galaxies. 
In semi-analytical calculations of the passage through a cuspy (p ~ density profile) 
primary potential, of a smooth secondary, also with a cuspy potential, one finds |6^ 



that, for reasonably elongated orbits, the tides of the primary pump in considerably 
more energy into the secondary than its own binding energy. Hence, secondaries such 
as galaxies dissolve at first passage. It remains to be seen whether these secondaries 
are completely disrupted or whether energy exchanges within it allow a dense core to 



survive. Recent simulations ||101|| of multi-particle secondaries moving along elongated 
orbits in a given isothermal potential show that tidal disruption is relatively slow, 
although it is not clear whether their "isothermal" primary is the singular p ~ 
model or the non-singular isothermal (that thus has an homogeneous core), for which 
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one indeed expects softer tidal effects. 

In any event, overmerging is much less evident in simulations where gas is included 
102| , |103|| , presumably because the gas sinks to the bottom of the galaxy potential 



wells and deepens these wells, which thus avoid merging with one another. The 
latter study emphasizes the importance of star formation in their gaseous galaxies in 
reducing overmerging. 

Dense groups of galaxies witness rapid merging and coalesce into a single elliptical 
galaxy [|10|, |T0|, H, |10|] . A detailed comparison of the results on groups ||107|| showed 



that the different numerical studies of groups produced comparable rates of merging. 
The differences arise in part from the initial conditions and from the fact that the 
merger rate in dense groups is decreased when most of the dark matter is distributed 



in a common envelope, rather than in individual galaxy halos ||108| , |105| , |90| , |109 
Indeed, with large individual halos, galaxy merging is direct, while the presence of 
an important intergalactic background causes galaxies to dissipate orbital energy by 
dynamical friction, suffer orbital decay and finally merge together at the group center 
[pO|] . It may be that overmerging is occurring in all these dense group simulations, 
and that with the inclusion of gas dynamics, the galaxies in dense groups may survive 
longer. 

5.2 Cosmological evolution 

Spherical top-hat cosmology (§ p.2|) , can be applied to a galaxy system bathing in an 
empty or uniform universe. An homogeneous isolated system should see its size evolve 
as shown in Figure la. It first follows the Hubble expansion, then decouples from 
this expansion and turns around, collapses and subsequently virializes. If a system 
is in dynamical equilibrium, one can apply the virial theorem, and so derive a virial 
mass Mvir = RV^/G. While the time evolution of the radius R, velocity dispersion 
V, ratio of virial mass to true mass Mvir/M and ratio of crossing time to running 
time tcr/t = R/iyt) parameters have been studied before ||110|| (though simple N- 



body simulations), and an analytical version is shown in Figures lb, Ic, and Id, one 
gains considerable insight in plotting the evolution of M^ir/M versus tcr/t as is done 
in Figure 2a. The dotted track is for groups made of point mass galaxies, while the 
solid track is for extended galaxies, which reach a terminal velocity at group collapse 
(because the smoothed potential is flat at the center), and after virialization, dissipate 
their orbital energy by dynamical friction against their common massive halo (merged 
from their individual halos after group collapse). 

The solid track in Figure 2a can be considered as a fundamental track that galaxy 
systems should follow, similar in concept to a Hertzprung-Russell diagram for cosmo- 
logical galaxy systems. Unfortunately, true masses of galaxy systems are not known. 
To compare with observed parameters, we must make an assumption on the true mass, 
and the simplest one is to assume that the true mass follows light, i.e., M/L = est. 
In Figures 2b, c, and d, we plot the observable parameters, M^-^/Lb vs. H^tcr, for 



loose groups |111| of different multiplicities, and superpose the theoretical evolution- 
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Figure 1: Time evolution of bias in observed radius (a), velocity dispersion (b), 
mass (c), and crossing time (d), relative to virial equilibrium (VE), where Tta is the 
turnaround time. The dotted curves show the evolution for point-masses, while the 
full curves show the effects of softened potentials and orbital energy dissipation by 
dynamical friction (starting at t — 3 Tta) ■ 



ary track, adjusting the y-axis with the high multiplicity groups of Figure 2b, while 
the scaling is imposed by theory. 

The high-multiplicity groups fit the theoretical tracks very well. A one proceeds 
to lower multiplicities, the statistical noise in the mass-to-light ratio and crossing time 
estimates increases, but so does the probability for chance alignments, which make 
the groups appear smaller while conserving on the average their velocity dispersion. 
Although precise assignments of group cosmo-dynamical states is difficult because 
of statistical noise, one can nevertheless get a handle on which groups are unbound 
(above theoretical track), which are still in their expansion phase (upper-right handle 
of track), which are near turnaround (lower-right handles of track), which are col- 
lapsing (central handle), which are near maximum collapse (first lower- left handle). 
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Figure 2: Mass assuming dynamical equilibrium, scaled to total mass (a) or total blue 
luminosity (b, c, and d), versus crossing time (in units of the age of the Universe for 
f2 = 1, while for Q = 0.2 the points should be displaced to the left by 0.1 decade). 
The polygons (b, c, and d) represent the loose groups. The thin curves are the 
theoretical point-mass evolutionary tracks, while the thick curves are the same for 
softened potentials and allowing for orbital energy dissipation after virialization. In 
(b, c, and d), these curves are scaled to mass-to-light ratios assuming that all groups 
have a true M/ L = 440 h. The typical error bar for each data point is shown. 



and those that are virialized (second lower- left handle). The theoretical fundamental 
track in observable space M^ij./L vs. tcr/^o thus represents a slice through the funda- 
mental surface (which is curved) of groups, where the third axis can be the scale of 
the system, such as its total optical luminosity. In contrast, an empirical fundamental 
plane has been reported [ |112[ for loose groups, similar to that of elliptical galaxies 
and clusters (see § p^ ). 

The true M/ L is obtained by extrapolating to the early virialized state (before dis- 
sipation of orbital energy, which occurs at nearly constant velocity dispersion since the 
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common halo should have near constant circular velocity). The loose groups ||111|| then 
have Mtrne/L = 440 h, much higher than the median M/L = 130 h, for the groups 
of > 4 members (the mass estimate used here is the median of the non-weighted 
virial, weighted virial, and projected masses). In other words, the mass-to-light ra- 
tios of groups are severely underestimated because most groups are still relatively near 
their turnaround phase. This points to f2 ~ 0.3 obtained by extrapolating the galaxy 
luminosity function [|113|| to (M/L)ciosure = 1560 /i. 

In any event, no groups in the loose group catalog has yet completed its collapse, 
not even the Virgo cluster included in the catalog, whose outer members are still 
collapsing onto the virialized core. The dynamical youth of loose groups suggests a 
high value (close to unity) for Qq. Details of this analysis are still in preparation. 

The fundamental tracks for groups can been computed in a hierarchical binary 
approximation, with a galaxy merging with an already merged pair. The precise 
position of the track depends slightly on the mass ratio of the merging pair. Any 
system can be represented at a given hierarchy. For example, the NGC 2300 group, 
in which X-rays from the intergalactic gas were first discovered ||114|| , is a loose triplet 
that has just turned around from its maximum Hubble expansion, while the tight 
binary around which the X-rays are roughly centered appears on the evolved part of 
the fundamental track (late collapse or rebound). The fundamental surface analysis 
can be easily applied to pairs, but the error bars in the position of a pair on the 
diagram are relatively important, and a detailed statistical analysis is required. Con- 
versely, clusters suffer little from statistical uncertainties of the observed parameters, 
but when a subcluster is merging into a cluster, the fundamental surface analysis on 
the primary cluster may break down as the secondary can exert non- negligible tidal 
stresses that alter the collapse of the primary. 

In the previrialization scenario, structures acquire angular momentum at maxi- 
mum expansion, and do not fully collapse to reach a rapid virial equilibrium. The 
fundamental track, should then be slightly altered, as it will follow the standard track 
at early epochs (and high observed crossing times, to the right of Fig. 2), and up to 
tcr/t ~ 0.03, then slightly rebound to the virialized horizontal wing (the thin hori- 
zontal track to the left of Fig. 2). If the previrialized system coalesces by dynamical 
friction, its track will then follow the soft collapse track, and its overall behavior 
will be similar to the standard track. However, a coalescence is somewhat akin to 
a collapse, but slower. Hence, the the decrease of the track to the lower left of the 
diagram is delayed with respect to the standard track. 

In Figure 3 are shown the locations of compact groups |12] in the fundamental 
track diagram. The groups on the upper left of the diagram are high velocity dis- 
persion groups close to the fundamental track. Conversely, the groups on the lower 
right of the diagram are low velocity dispersion groups too far from the fundamental 
track to be explained by standard evolution. Instead, they are best explained as 
arising from chance alignments of galaxies along the line of sight, either within larger 
loose groups ||115|| that are still collapsing, or within long (~ 7 Mpc) filaments 



REFERENCES 



21 



|p.l6|| , often seen in cosmological simulations, and some of which may be expanding 
(the spherical equivalent being the upper right track). Thus, the higher velocity 
dispersion half of this sample of compact groups are real 3-D dense systems, either 
terminating their collapse or rebounding from it, or in their final coalescence stage, 
while the lower velocity dispersion half are caused by chance alignments. If previri- 
alization is not followed by rapid coalescence, then as much as 90% of the compact 
groups would lie too far from the fundamental track to be real and would be ascribed 
to chance alignments. This is inconsistent with the detection of X rays, probing the 
potential of truly dense groups, in roughly half of these compact groups [|1171| (typ- 
ically the high velocity dispersion ones ||117| , |118|| ). If coalescence does occur after 
previrialization, the data of Figure 3 is not good enough to distinguish between the 
two scenarios, but a richer compact group catalog, or the use of X-rays to infer in a 
more secure fashion the velocity group dispersion should help settle between them. 



1000 




0.01 



0.1 

t../t, 



Figure 3: Same as Figures 2bcd for compact groups {crosses] 
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